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Abstract 



1 State of Art 



CREME96 and GEANT4 are two well known particle 
transport codes through matter in space science. We 
present a comparison between the proton fluxes out- 
going from an aluminium target, obtained by using 
both tools. The primary proton flux is obtained by 
CREME96 only, covering an energy range from MeV 
to hundreds GeV with the same result in both cases. 
We studied different thickness targets and two different 
GEANT4 physics lists in order to show how the spectra 
of the outgoing proton fluxes are modified. Our findings 
show consistent agreement of simulation data by each 
tool, with regards both GEANT4 physics lists and every 
thickness target analysed. 



Introduction 

The scientific world makes use of several simulation 
codes to study the interactions between radiation and 
matter, as can be see in the overview reported in sec. [T] 
In our study, we focused on two of these, CREME96 
and GEANT4 both described in sec. [2] These two toolk- 
its are characterized by different physics, both illus- 
trated in sec. [5] In order to realize this study we imple- 
mented a specific simulation set-up, as defined in sec.|4j 
and we used a statistic data analysis explained in sec. [5] 
Our results, reported in sec. [6l show the agreement be- 
tween CREME96 and GEANT4. As we suggest in El 
these promising results have encouraged us to under- 
take further comparisons using different types of pri- 
mary particles. 



Nowadays there are several particle transport codes 
through matter, for example CREME86, CREME96Q, 
FLUKA0, GEANT30, GEANT40, MARS0, 
MCNP01, MCNRXEl, PHITS. 

These tools use different languages, but most of 
the codes are designed in fortran; f77 for GEANT3, 
FLUKA, PHITS; f90 in the case of MCNPX and f95, 
for MARS. The only one of these examples in C++ is 
GEANT4. Source codes for CREME86 and 96 are not 
freely available. 

Each is developed by a different community. The 
owner of CREME, 86 and 96, is the USA Naval Re- 
search Laboratory. On the other hand MCNPX is de- 
veloped by LANL. The GEANT4 community includes 
CERN, ESA, IN2P3, PPARC, INFN, LIP, KEK, SLAC, 
TRIUMF. FLUKA support is made by CERN and INFN 
researchers and MARS is a FNAL laboratory produc- 
tion. Finally, PHITS code is managed by JAEA, RIST, 
GSI and Chalmers University. 

Because of these differences, working simultane- 
ously with all these codes is a very hard task. Therefore 
a scientist usually carries out research on one particular 
code and sometimes makes a comparison between their 
own software and another code. Examples of this kind 
of comparative study [ 8 1 have already beeen undertaken 
for FLUKA and GEANT4 and [9| for CREME 86 and 
GEANT4. 
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2 Simulation tools considered 

As previously stated (sec. [TJ, there are many tools sim- 
ulating the passage of particles through matter. In our 
study we focused on two of these: CREME96 and 
GEANT4. 

• Cosmic Ray Effects on Micro-Electronics 
(CREME) is a toolkit developed to study the 
effects of an ionizing radiation environment on 
the electronics inside a satellite, i.e. Single Event 
Upsets (SEU) due to solar heavy ions iflOl . 

The CREME96 user manages the simulation by 
point and click method, by means of an easy-to- 
use graphical interface avaiable on-line at Q]. 

CREME96's principal features, for instance the 
physical models used and their limitations, are re- 
ported in ifTTI . These models and their constraints 
are an important element of our study, for this rea- 
son we will describe them in more detail below. 

• GEANT4, GEometry ANd Tracking (particles), is 
a Monte Carlo toolkit for simulating the passage of 
the particles through matter, avaiable on-line at [4 ] . 

GEANT4's first production release was completed 
in December 1998 fl2l . Unlike its ancestor 
GEANT3, in fortran 77, GEANT4 is designed in 
C++. 

In GEANT4 there is not a default physics list, 
therefore the user must define his/her specific 
physics. The features of our GEANT4 physics list 
are explained in sec. [3] 

Moreover, in GEANT4 there is no pre- 
implemented definition of particle flux. In 
order to overcome such a problem, QinetiQ 
developed MUlti-LAyered Shielding Simulation 
Software (MULASSIS [CD), a GEANT4 inter- 
facing tool, which also contains other features to 
work with this Monte Carlo code. 

We did not consider MULASSIS for two reasons. 
Firstly, because we prefered to make a direct com- 
parison between CREME96 and GEANT4 with- 
out being compelled to also take into account the 
additional aspect of an intermediate programme. 
Secondly, we needed to input both CREME96 and 
GEANT4 radiation transport calculation packages 
with exactly the same particle flux output (the 
same particle flux output produced as a text file 



by CREME96). Conducting the same study us- 
ing MULASSIS would have added extra compli- 
cations. 

For these two reasons we did not use MULASSIS, 
deciding instead to work with a simple code suit- 
able for our research purposes. 

3 Physics in the two simulators 

CREME96 and GEANT4 were developed by two dif- 
ferent scientific communities in different time periods. 
As a result their physics show a number of differences. 

• The physics of CREME96 includes models of the 
Galactic Cosmic Rays (GCRs), Anomalous Cos- 
mic Ray (ACRs), Solar Particle Events (SPEs) and 
a geomagnetic transmission calculations. These 
models are used for building a flux of particles 
about an orbit defined by the user. The proton flux 
used in our study is generated by CREME96 only, 
because it is not possible for GEANT4 to do it. 

Another aspect of CREME96's physics model is 
the nuclear transport physics through the shielding 
material. This tool takes into account both ioniza- 
tion energy loss (in the continuous slowing down 
approximation) and nuclear fragmentation. It im- 
plements stopping power and range-energy rou- 
tines and uses semi-empirical energy-dependent 
nuclear fragmentation cross-sections lfTTI . 

The physics in CREME96 has three limitations. 
First of all, it does not consider neutron produc- 
tion in the shield, because these particles do not 
cause damage to electronics of spacescraft. Sec- 
ondly, it does not include target fragments. Finally, 
aluminium is the only type of shielding material. 

• CREME96 can describe both the external en- 
vironment (no shield) and the internal satellite 
(shielded) via flux of particles with Z range from 
1 to 28, without requiring its user to define or to 
implement their flux definition. 

• With regards GEANT4, in this study we used two 
types of physic list, named type A and type B, in 
order to see the differences between the two sim- 
ulation codes, but also between our two particular 
choices. 
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As far as the electromagnetic physics list was con- 
cerned both types included the following mod- 
els: photon-epdl; electron-eedl; positron-standard; 
ion-ICRU; muon-standard and decay. 

The difference between the two types is in the 
hadronic physics list selected, Binary for type A 
and Bertini for type B. 

The physics implementation used in our study 
depended on many elements; in the case of 
GEANT4, for example, its patch and its library 
version. In our particular study we used Geant4 
version 8.1 patch 1 and CLHEP library version 
1.9.2.3.. 

• As previously stated, GEANT4 does not have a de- 
fault definition of particle flux, therefore we had to 
solved it in the first instance. In our implemen- 
tation of the particle flux (4>) in order to make a 
comparison for CREME96 <t> was defined as fol- 
lows IT4ll : 

Def.: 1 (Flux) "amount of radiation crossing a 
surface per unit of time, often expressed in "inte- 
gral form" as particles per unit area per unit time 
(e.g. electrons cmT 2 s~ x ) above a certain threshold 
energy. 

NOTE The directional flux is the differential 
with respect to solid angle (e.g. particles 
cm~ 2 steradian^ 1 s~ l ) while the "differential" 
flux is differential with respect to energy (e.g. par- 
ticles cm~ 2 MeV^ 1 S~ ! ). In some cases fluxes are 
also treated as a differential with respect to Linear 
Energy Transfer (. . . )." 



4 Simulation set-up 

The considered geometric configuration and the mate- 
rial set up of the experiment were reproduced in both 
simulators in the same way and consisted of two ele- 
ments: the beam and the target. 

• The Beam. A punctiform beam made of protons 
was aimed at the face of the target, towards its cen- 
ter. 

As previously stated (ref. sec. [3]) the number of pri- 
mary particles generated during the simulation can 
not be explained in this study, because the flux <t> 
(ref. def. [U was normalized by primary particle 
number. 

The energy weight of a single primary proton de- 
pends on the flux used with input for both simula- 
tions. This flux is generated only by CREME96, 
as explained in sec. [3] For the purpose of our pa- 
per, it derived from the choice of a typical ISS or- 
bit thus defining its component environmental fea- 
tures (see fig. IB- 

• The target. This was made of aluminium because, 
as previously stated (ref. sec. 0), CREME96 can 
only functions with this material. 

The target thickness is made to vary along the line 
normal to the beam direction, from 1 -5 cms at 1 cm 
intervals, in order to observe how the outgoing 
proton flux was modified. These thickness values 
were chosen because they are typically used when 
modelling the satellite shielding action in a simple 
way, e.g.: FOTON-M 3 [15] can be approximated 
with a 3 cm layer of equivalent aluminium. 



In this paper, according to def[U the measurement 
unit of the flux used in the input file is: 



particles 



m 2 ■ s ■ sr- MeV jnucl 



(1) 



where m is meter, s is second and sr is steradian 
and it is normalized by the number of primary par- 
ticles generated during the simulation. 



5 Simulation Analysis 

We limited our study to the proton fluxes both for simu- 
lator input and simulator output. In particular, we clar- 
ified that in the case of output in the particle flux there 
were both primary protons going out of the target and 
secondary protons produced in the shield. 
CREME96's fluxes were not normalized by the number 
of primary particles, therefore in order to make a com- 
parison we had to first perform normalization. For this 
tool we calculated the primary number by integration 
of the proton flux used for simulation input. The result 
was 8.894 ■ 10 4 particles, and is reported in fig. 03 
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The data output of the two simulators did not have the 
same statistics. 

CREME96 data did not have the graphical error bar, 
whereas GEANT4 simulations received the following 
statistic: 

x-axis error bars were constant and are equal to the 
binning chosen by the user during the implementation 
code. In this study we have defined a 1 MeV/Nucleon 
bin. 

The y-axis error bars were variable. The dimension 
of this bar type followed the Poisson statistic i.e.: 



where N is the number of entries w.r.t. the relative 
energetic bin. 

6 Results and comparison 

Our results are plotted in fig|2lfTT1below shielded proton 
fluxes are reported using two kinds of scale, decimal 
(figs. El H |6] IdOl) and logaritmic (figs. EE] E3IID- 
In the first case we can clearly observe the peak corre- 
sponding to particle energy with maximum transit prob- 
ability. 

In the logarithmic version it is possible to view the 
energetic range of the particles and how the flux de- 
creases. In this representation the error bars of the en- 
ergetic particles are very wide because there are very 
low entries for the higher energies. We maintain that 
the number of primary particles needed to obtain good 
error bars for high energy particles (w.r.t. our primary 
flux) should be more than 10 10 . This type of simulation 
is too expensive for our computational infrastructures, 
therefore we had to limit our simulation to 10 6 primary 
particles. 

For all shielding thicknesses, from 1 cm to 5 cm, we 
can observe a good agreement between CREME96, and 
GEANT4 simulation data, for both physics list type. 
Sometime the overlapping of the different data is so 
good that we can not see the CREME96 data. In fig- 
ures where a decimal scale is used, we can distinguish 
the visible difference on the peak. This effect is due 
to different nature of the simulation code output data: 
CREME96 reports a continuous curve, GEANT4 shows 
the information data in a histogram therefore the flux 
shape is not a continuous curve. 



7 Conclusion and future work 

We made a comparison between CREME96 and 
GEANT4, two particle transport codes through mat- 
ter, in order to measure an aluminium target's outgo- 
ing proton fluxes and up until now evaluated using only 
CREME96. To test GEANT4, we worked activating 
two different kinds of physics lists. Our results show 
a good agreement between the output of the two simu- 
lators for each GEANT4 physics list considered. 

These positive results encourage us to carry out more 
research in the same field in the future. A possible evo- 
lution of this work is, in fact, the analysis using other 
primary particles, since CREME96 can consider parti- 
cles with Z varying from 1 to 28. 
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Figure 1 : Primary proton fluxes, plotted in logarithmic 
scale. 
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Figure 2: Proton fluxes behind 1 cm of aluminium 
shielding, plotted in decimal scale. 




Figure 3: Proton fluxes behind 1 cm of aluminium 
shielding, plotted in logarithmic scale. 
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Figure 4: Proton fluxes behind 2 cm of aluminium 
shielding, plotted in decimal scale. 
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Figure 5: Proton fluxes behind 2 cm of aluminium 
shielding, plotted in logarithmic scale. 
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Figure 6: Proton fluxes behind 3 cm of aluminium 
shielding, plotted in decimal scale. 
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Figure 7: Proton fluxes behind 3 cm of aluminium 
shielding, plotted in logarithmic scale. 
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Figure 8: Proton fluxes behind 4 cm of aluminium 
shielding, plotted in decimal scale. 
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Figure 9: Proton fluxes behind 4 cm of aluminium 
shielding, plotted in logarithmic scale. 
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Figure 10: Proton fluxes behind 5 cm of aluminium 
shielding, plotted in decimal scale. 




Figure 11: Proton fluxes behind 5 cm of aluminium 
shielding, plotted in logarithmic scale. 
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